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ABSTRACT 

We present results from an investigation of the dynamical behavior of 
buoyant magnetic flux rings in the radiative interior of a uniformly rotating 
early-type star. Our physical model describes a thin, axisymmetric, toroidal 
flux tube that is released from the outer boundary of the convective core, and is 
acted upon by buoyant, centrifugal, Coriolis, magnetic tension, and aerodynamic 
drag forces. We find that rings emitted in the equatorial plane can attain a 
stationary equilibrium state that is stable with respect to small displacements 
in radius, but is unstable when perturbed in the meridional direction. Rings 
emitted at other latitudes travel toward the surface along trajectories that 
largely parallel the rotation axis of the star. Over much of the ascent, the 
instantaneous rise speed is determined by the rate of heating by the absorption 
of radiation that diffuses into the tube from the external medium. Since the time 
scale for this heating varies like the square of the tube cross-sectional radius, for 
the same field strength, thin rings rise more rapidly than do thick rings. For a 
reasonable range of assumed ring sizes and field strengths, our results suggest 
that buoyancy is a viable mechanism for bringing magnetic flux from the core 
to the surface, being capable of accomplishing this transport in a time that is 
generally much less than the stellar main sequence lifetime. 



Subject headings: Stars: magnetic fields Stars: early-type — Stars: interiors — 
Stars: radiative envelope — Stars: surface fields — 



^High Altitude Observatory, National Center for Atmospheric Research, P.O. Box 3000, Boulder, CO 
80307-3000; kmac@hao.ucar.edu 

^Department of Astronomy, University of Wisconsin, 475 N. Charter St., Madison WI 53706-1582; 
cassinelli@astro . wise, edu 



- 2 - 



1. Introduction 

Although magnetic fields and related activity are believed to be nearly ubiquitous 
among stars of solar and lower mass, it is less certain that magnetism is a common 
characteristic of stars more massive than the Sun. To date, definitive detections of fields 
in stars with masses ^ 1.5 Mq have, for the most part, been made for objects having 
anomalous chemical abundances (e.g., the chemically peculiar A and B stars; Landstreet 
1992; Donati 1998). Recently, however, observations of cyclic variabihty in the properties 
of winds from luminous OB stars have been interpreted as evidence for the presence of 
large-scale magnetic fields in the surface layers and atmospheres of these objects (see, e.g., 
Kaper & Henrichs 1994; Kaper et al. 1997; Prinja 1998). These inferences have been 
bolstered by the unambiguous measurement of a weak (~ 360 G) field in the chemically 
normal Bl Ille star /3 Cephei (Donati et al. 2001). These results suggest that magnetic 
fields of moderate strength might be more prevalent among hot stars than had previously 
been thought (see Charbonneau & MacGregor 2001, hereafter Paper I, for a more complete 
discussion) . 

At the present time, the origin of magnetism in massive stars is not well understood. 
It is unclear whether the inferred fields are fossil in nature, or are instead the product of 
dynamo activity in the stellar interior (see, e.g., Parker 1979). The latter explanation, if 
appropriate, raises additional interesting questions concerning the site of magnetic field 
generation inside hot stars. For stars having spectral types O and B, convection occurs in 
the innermost portion of the core, not in the outer envelope as in low-mass stars. Because 
convection is thought to be necessary in order that field regeneration via the so-called 
a-effect take place, it follows that a hot star dynamo should be located deep within 
the stellar interior. Within the context of mean-field electrodynamics, the properties of 
kinematic core dynamo models for early-type stars have been investigated by Levy & Rose 
(1974) and Schiissler & Pahler (1978), and, more recently, in Paper I. 

If the magnetic field of a hot star is produced by dynamo action in the convective 
core, then a mechanism for transporting the field to the stellar surface must be identified. 
The finite electrical conductivity of the envelope leads to the outward diffusion of any 
fields contained therein, but only over an extended period of time. Estimates indicate 
that for stars more massive than a few solar masses, the resistive diffusion time across 
the radiative interior exceeds the main sequence lifetime tMS (Schiissler & Pahler 1978). 
Another possibility is that dynamo fields are advected from the core to the surface by 
rotation-induced meridional circulation (see Paper I). For a star of mass and equatorial 
rotation speed Vroti a rough estimate for the circulation time tf. as a fraction of Ims is 
itc/tus) ^ 9 X 10^ (M,/Mq)°-s (^^ro^/l km s"^)-^ (Kippenhahn & Weigert 1994). In the case 
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of a star with = 10 Mq and Vrot — 150 km s"-*^, (tc/tMs) ~ 0.25, indicating that even for 
relatively rapidly rotating stars, a time t\is is required to bring the field produced by the 
dynamo to the surface by this mechanism. In addition, the results of Paper I suggest that 
the circulatory flow generated by extremely rapid rotation is capable of interfering with the 
operation of a core dynamo. 

Alternatively, in the Sun, magnetic flux emerges at the photosphere in the form of 
flbrils or flux tubes. The fleld is thought to assume this form in or near the dynamo domain 
at the bottom of the convective zone. The subsequent rise of a flux tube to the solar 
surface is driven by buoyancy, a consequence of the reduced density inside a tube that is in 
mechanical and thermal equilibrium with the surrounding, adiabatically stratifled, fleld-free 
gas. If the formation of flux tubes is a process that is not specific to fields in the Sun, then 
a dynamo-generated field at the base of the extended radiative envelope of a hot star might 
also develop fibril structure. In this case, the buoyant force might likewise enable fiux from 
deep in the interior to reach the surface in a time that is shorter than evolutionary time 
scales. An important distinction, however, is that unlike the the solar convection zone, the 
stratification within the radiative interior of a hot star is sub-adiabatic. The motion of 
buoyant magnetic elements in a thermodynamic environment like that in the envelope of an 
upper main sequence star has been studied by Gurm & Wentzel (1967) and Moss (1989). 

In the present paper, we adopt the premise that hydromagnetic dynamo activity takes 
place inside the convective core of a rotating massive star, in a manner similar to that 
described in Paper I. We furthermore presume that the flelds so-produced naturally form 
into discrete, toroidal flux tubes that remain in total pressure equilibrium with the external, 
unmagnetized stellar interior at all times. Using this conceptual framework in concert with 
the physical model described in §2, we determine the time-dependent position, velocity, and 
thermodynamic properties of a buoyant flux ring that begins its outward motion from a 
specified location on the core boundary. Results pertaining to the dynamics and rise times 
of rings with a variety of initial field strengths and cross-sectional radii are presented in §3. 
Our findings regarding the dynamical behavior of magnetic rings in a radiative environment 
and the efficacy of buoyancy as a fiux transport mechanism in hot stars are summarized in 
§4. 

2. Model 

To investigate the buoyant transport of magnetic flux through the radiative interior of 
a hot star, we consider a spherical star of mass and radius that contains a central, 
convective core of mass Mc (< M*) and radius Rc (< i?*). The star is assumed to rotate 
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uniformly with angular velocity Q around an axis that coincides with the polar axis of a 
spherical coordinate system (r,6,(f)). In the radiative region external to the core, wc study 
the dynamics of an isolated concentration of azimuthally directed magnetic field (i.e., a 
flux tube) that at any given time takes the form of a circular ring, symmetric about the 
stellar rotation axis. A similar configuration has been utilized by Choudhuri & Gilman 
(1987) to examine the influence of rotation on the rise of flux tubes in the solar convection 
zone. We assume that the cross-sectional profile of the tube is a circle whose radius a is 
much smaller than either the radius of curvature of the ring or the local scale height of the 
surrounding, unmagnetized, radiative envelope. The assumed geometry for the rising fiux 
tube is illustrated in Figure 1. These approximations are here used to study the dynamics 
for the rise of magnetic fiux to the stellar surface, and for comparisons of rise times with 
the main sequence lifetime of the star. 

As it is described above, the magnetic ring satisfies the criteria for applicability of the 
thin fiux tube approximation (e.g., Defouw 1976; Spruit 1981), and its internal properties 
can be reasonably taken to be uniform throughout. The ring is assumed to contain material 
with mass density p, so that its total mass is m = pV, where V = 27rr sin 6 ■ vra^ is the ring 
volume. It follows from conservation of mass that if a ring with cross-sectional radius Oq is 
released from position (ro, 6*0) at time t = 0, its density at some t (> 0) is 

/ao\ V'^o sin 6'o\ 
\ a J \ r sm B J 

where a and (r, 9) are the tube radius and position at the later time. Likewise, the fiux ring 
is presumed to be untwisted and to contain a toroidal magnetic field of the form B = B 
with B a constant. Conservation of magnetic fiux then implies that B is related to its 
initial value Bq through the relation 

B = bJ"-^)\ (2) 



a 

A further consequence of the thin fiux tube approximation is that because the time required 
for a fast magnetosonic wave to transit the tube cross section is short, the total pressure 
inside the ring can be assumed to instantaneously equilibrate with that of the surrounding 
stellar material. The quantitative expression of this lateral mechanical equilibrium condition 
is 

B^ 

Pe^P+:^, (3) 

where p — pkT/ p, is the gas pressure and the subscript e refers to conditions in the external, 
field-free medium on the periphery of the tube. Note that by using the ideal gas law for 
p and Pe, we are neglecting the radiative component of the total pressure. This omission 
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restricts our analysis to stars less massive than about 10 Mq, for which radiation pressure 
is a minor contributor to the support of the deep stellar interior. In addition, we assume 
that the mean molecular weight ^ has the same constant value throughout the star. Our 
model thus describes conditions within a young, chemically homogeneous, main sequence 
star; it does not apply at later evolutionary stages when a flux tube that originates in the 
vicinity of the core-envelope interface may contain chemically processed material. 

In a frame of reference that rotates with the stellar angular velocity, the uniform flux 
ring moves coherently in response to the axisymmetric forces applied to it; that is, each 
mass element within it has the same velocity 

dr d9 . /, d0 

U^Ur Br + Ug eg + Ud> = — Br + r — Be + r Sm — 60. (4) 

at at at 

The velocity components {ur, ug, u^) and ring position (r, 9) can be obtained as functions of 
time by integration of the ring equation of motion. 



dt \Pe + P, 



GM{r) 



Br — fl^ r (sin^ 9 Br + sin 9 cos 9 Bg) 



(e, + cot 9 Be) - Cn ( (5) 



Attt (pe + p) \Pe + pj Tra 

where M(r) is the internal mass distribution of the star, Cd is the drag coefficient, and 
u± = UrBr + uq Bq is the trausvcrsc velocity. The individual components of the inertial and 
Coriolis terms on the left-hand side of equation (5) have the explicit forms 

\ — — h 2 iZ X u = — 2il U(k sm 9, (6a) 

\dt J ^ dt r ' y J 

(^ + 2nxu] ^^ + ^-^cot9-2nu,cos9, (6b) 
\dt J Q dt r r 

(^ + 2f2xul =^ + ^ + ^cot^ + 2Q (ii,sin^ + ii,cos^). (6c) 
\ dt }^ dt r r 



The equation of motion (5) and its many variations have been used to treat the 
dynamics of buoyant flux tubes in the Sun's convective envelope by numerous authors, 
including Choudhuri & Oilman (1987), Ohoudhuri (1989), Moreno-Insertis, Schiissler, & 
Ferriz-Mas (1992), Oheng (1992), Fan, Fisher, & DeLuca (1993), Ferriz-Mas & Schiissler 
(1993, 1994), Fan, Fisher, & McOlymont (1994), Oaligari, Moreno-Insertis, & Schiissler 
(1995), and Fan & Fisher (1996), among others. The three terms on the right-hand side of 
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equation (5) describe, respectively, the accelerations produced by: (i) the buoyant force, 
including the effect of the centriftigal reduction of the local gravitational acceleration: (ii) 
the magnetic tension force, arising in response to outward, stretching motion of the ring in 
the plane perpendicular to the rotation axis of the star; and (in), the aerodynamic drag 
force, derived from that exerted on a straight, circular cylinder immersed in a steady flow 
with upstream velocity directed perpendicular to the axis of symmetry (Goldstein 1938). 
The appearance of the factor (pe + p)~^ rather than in the expressions for each of the 
accelerations reflects the fact that we have made provision for the so-called virtual (or 
enhanced) inertia of the ring, the hydrodynamical resistance to acceleration experienced by 
a rigid body that moves through a fluid (see, e.g., Batchelor 1967). Considerable discussion 
has been devoted to the reasons for including (or not including) this effect in the equation 
describing the balance of forces acting on a flux tube (see, e.g., Moreno-lnsertis, Schiissler, 
& Ferriz-Mas, 1996; Fan & Fisher 1996, and references therein). As will become clear in 
subsequent sections, the overall picture of flux ring dynamics that emerges from solutions of 
equation (5) depends little on the presence or absence of virtual inertia in any of the terms. 

The difference between the nature of the energy balance that prevails within the 
envelope of a massive star and that within the convection zone of the Sun leads to signiflcant 
differences between the dynamics of flux tubes located in the two regions. The outer 
portion of the solar interior is adiabatically stratified, a consequence of efficient convective 
energy transport therein. Because of this, an adiabatic fiux tube that is initially buoyant 
will remain so during the course of its rise toward the surface. Alternatively, the envelope 
of a hot star is sub-adiabatically stratified, a consequence of the radiative equilibrium 
conditions that exist throughout it. As a result, the density deficit of an initially buoyant, 
adiabatic tube will diminish as it rises, causing the upward buoyant acceleration to decrease 
as well. When adiabatic cooling causes the temperature contrast of the tube relative to its 
surroundings to attain the value (Tg — T)/Te = (1 + (3)~^ where (3 = {STrp/B"^), p = p^, and 
the buoyancy of the tube is reduced to zero. 

The discussion of the preceding paragraph indicates that within the envelope of 
a massive star, the non-adiabatic, radiative interaction between a flux ring and its 
environment will play an important role in determining the magnitude of the buoyant 
acceleration the ring experiences, and thus, the time required for it to emerge at the surface 
of the star (see, e.g., Parker 1975; Moreno-lnsertis 1983; Fan & Fisher 1996). In the present 
model, the heating associated with the diffusion of radiation from the external medium into 
the tube is described by the thermal energy equation (see Appendix A) , 

dS _ 32 (7 - 1) a AT 

dt 3 Ke Pe P 0? Te 



(7) 
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where AT = {Tg — T) , and 5" is the entropy-like quantity 



5 = In 




(8) 



In equations (7) and (8), 7 is the ratio of specific heats, a is the Stefan-Boltzmann constant, 
and Ke is the Rosseland mean opacity in the material surrounding the tube. In the absence 
of radiative heating (i.e., when q — 0), equation (7) indicates that S is constant, and the 
tube behaves adiabatically. 

A general question concerning the flux tube approach is how far we can consider the 
flux tube to be thin in the thermal sense. We assume the heating that arises from the 
inward diffusion of radiation affects the tube interior uniformly. In reality, the fact that the 
heating time scale and the tube radius are both flnite implies a differential response by the 
tube to the heat input. In particular, if the time scale for the dynamical adjustment of the 
tube (~ the sound crossing time) is much shorter than the diffusion heating time scale, then 
the temperature in the outermost portion of the tube can become elevated relative to the 
center, and the resultant enhancement of the buoyant acceleration can cause a deformation 
of the tube. Unless checked by some additional conflning influence (e.g., twisted flux tube 
fields), this tendency might ultimately lead to the destruction of the fiux concentration. 
While quantitative treatment of this effect is beyond the scope of the present exploratory 
model calculations, its potential occurence underscores the desirability of performing 2D 
MHD simulations of fiux tubes in a radiative enviroment. 

We use the model consisting of equations (l)-(8) to determine the motion of an initially 
buoyant fiux ring in the radiative envelope of a hot star. The distributions of pressure (pe), 
density (pe), and temperature (Tg) are taken from a spherical, non-rotating stellar interior 
model for M* = 9 M©, kindly calculated for us by S. Jackson. As described in Appendix 
B, in order to facilitate numerical computations, we have derived an approximate, analytic 
representation of this model by assuming that pe, Pe, and Tg are related polytropically. For 
the purpose of these simulations, the star is assumed to rotate rigidly with a prescribed 
angular velocity Q that is sufficiently slow that departures from sphericity can be neglected. 

At time t = 0, a toroidal flux ring with specifled values of the cross-sectional radius qq 
and plasma beta f3o (= Sttpq/Bq) is released from the outer boundary of the convective core 
(ro = Rc) at co-latitude 9q (latitude Aq). We suppose that the ring is initially in thermal 
equilibrium with the ambient stellar material, Tq = T^q. Application of the mechanical 
equilibrium condition (3) then yields po = Peo Po/ (1 + Po) and po = Peo Po/ (1 + Po), so that 
(Peo — Po) > and the ring is buoyant. The initial field strength Bq and Alfven speed uao 
in the ring are given in terms of (3o and the external properties at the starting position 
according to Bq = [87rpeo/(l + Po)]^^"^ and uao = [2peo/(/3oPeo)]^'^^- 
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The velocity components {ur,uo,u^) are obtained as functions of t by numerical 

integration of the three components of equation (5); the position {r,9) of the axisymmetric 
ring follows by integrating the relations dr/dt = Ur, d9/dt — uq/t given in equation (4). 
In addition, we simultaneously integrate the thermal energy equation (7) to obtain S, 
which enables us to update the tube properties at each time step. In particular, from the 
definition given in equation (8), the pressure p within the tube can be expressed as 



p=p„e»(^y. (9) 



Similarly, if equation (1) is rewritten in the form 

a \ ^ ( Po\ ( tq sm. 6q 



Mo J \ p J \ r sin 9 
then substitution in equation (2) yields 



(10) 



b = bJ^)(iJ!^]. (11) 

\PoJ \ro sm 9oJ 

This result, together with equation (9), permits the pressure equilibrium condition (3) to 
be recast as 

,.(iL)\l(^)ViLy_(i,i)(M.o. (12) 

\PoJ Po\rosm9oJ \po J \ Po J \Peo J 

which is solved to obtain p/po, given the values of 5* and Pe that prevail at a particular time 
t and location {r,9). Equations (9), (10), and (11) are then used to derive p, a, and B, and 
T is determined via the relation T/Tq = (p/po)(Po/p)- 



3. Results 

We have used the thin flux tube physical model and solution method outlined in 
the preceding section to investigate the dynamical behavior of magnetic flux rings with 
10^ < Po < 10^ and 10"^ < (ao//ieo) < 10"^ where heo (= 2.65 x 10^° cm; sec Appendix 
B) is the pressure scale height in the stellar interior at the starting radius Tq = Rc- The 
maximum and minimum values adopted for the quantity /3o correspond to toroidal fields 
with strengths Bq — 1.82 x 10^ G and 5.77 x 10^ G, respectively. In the computational 
results described below, the parameters 7, Co, and ft were assigned the values, 7 = 5/3, 
Cd — 1, and il — 5.88 x 10~^ s~^, the latter value representing a surface equatorial 
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rotation speed of about 150 km s~^. We have examined the motion of buoyant flux rings 
following release from the core boundary at initial latitudes in the range 0° < Aq < 45° 
(i.e., 7r/4 < ^ Because the dynamics of rings in the equatorial plane (Aq = 0°) 

differ from those of rings with Aq > 0°, we treat the two cases separately in the ensuing 
discussion. 



3.1. Flux Rings in the Equatorial Plane 

Inspection of the equation of motion (5) reveals that the meridional component of 
each force acting on a flux ring that is initially at rest (u — 0) in the equatorial plane 
{6 = 7t/2) vanishes identically. The subsequent motion of such a ring therefore remains in 
this plane, and is governed by just the r- and 0-components of equation (5). Some of the 
properties of configurations of this kind are illustrated in Figures 2 and 3, which contain 
results pertaining to rings with — 10^. The pairs of panels that make up the three rows 
of Figure 2 show, respectively, the time evolution of: (i) the radial position of the ring, 
(Ar/ro) — {r — ro)/ro; (ii) the velocity components {ur,uo,U(j)) in units of the Alfven speed 
uao {— 7.49 X 10^ cm s~^) in the ring a,t t — 0; and (iii), the accelerations produced by 
each of the forces acting on the ring, in units of u\q/R^. In each panel, time is measured in 
units of the Alfven time i^, defined as Ia = R*/uao {— 3.42 x 10^ s ~ 4 days). 

The dynamical evolution of a flux ring with initial cross-sectional radius qq — 10~^ heo 
is shown in panels (A), (C), and (E) of Figure 2. Note that in this case, the outward 
expansion of the ring ceases after t !v Ia, its radial position remaining flxed at later 
times. In the course of evolving toward this apparent equilibrium configuration, the ring 
experiences a brief, initial period of accelerated motion in the +e,.-direction, followed by 
radial deceleration to a state of rest in which Ur = but 7^ 0. Examination of the 
ring force balance indicates that the initial expansion is driven (as expected) by buoyancy, 
while being opposed by the slightly smaller, inward-directed, magnetic tension force. The 
aerodynamic drag force, although small in magnitude, contributes to the damping of the 
oscillatory motions that are induced by the interplay between the larger buoyant and 
magnetic forces. 

The flux ring is assumed to corotate with the stellar interior at i = 0, so that its initial 
azimuthal velocity is li^o = in the reference frame that rotates with angular velocity 
n. For t > 0, the azimuthal motion of the ring is such that the angular momentum of 
the material contained within it is conserved. The buoyancy-driven increase in the radial 
position of the ring is therefore accompanied by a decrease in the rotation rate of the ring 
material. Hence, in the frame of reference that corotates with the star, the azimuthal 
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velocity of the ring is in the — e<^-direction, as can be seen in panel (C). According to 

equation (6a). this azimuthal motion produces an inward-directed Coriohs force that grows 
with increasing Ar until its magnitude, combined with that of the magnetic tension force, 
is sufficient to balance the buoyant force and prevent further radial expansion of the ring. 
A similar behavior has been observed in connection with models of equatorial magnetic 
flux rings located in the radiative layers just below the base of the solar convection zone 
(Moreno-Insertis, Schiissler, and Ferriz-Mas 1992; Ferriz-Mas 1996). 

In panels (A) and (C) of Figure 3, wc show how the thermodynamic state of a ring 
with (3q = 10^ and Oq = 10~^ heo evolves over time. Panel (A) shows the normalized 
differences between the pressure, density, and temperature in the tube and in the material 
surrounding it (i.e., Ap/p^, Ap/p^, and AT/Tg, respectively, where for any quantity Q, 
AQ = Qe — Q) as functions of time. Beginning from a state in which Tq — Tgo, the ring 
cools slightly relative to its surroundings during the earliest and most rapid portion of its 
limited buoyant ascent. As its outward progress stalls, radiative heating acts to restore the 
condition of thermal equilibrium between the ring and its external environment. In panel 
(C), we display the time evolution of the quantities S and q (see equations [7] and [8]), 
along with the corresponding histories of the instantaneous heating and rise times, defined 
as tn = (see Appendix A) and = \he/ur\, respectively. Note that S increases as the 
tube is radiatively heated, and approaches a constant value as — >^ 0, T — > Tg, and both 
tn and tji become large. 

To further explore the nature of the equilibrium seen in Figures 2 and 3 for fiux rings 
in the equatorial plane, we write the 0-component of the equation of motion (5) in the form 

r— + UrUs + 2urO,r = — (ru^ + fir^) = 0, (13) 
at at ^ 

where we have assumed ue = for 9 = 7i/2, and have used the fact that Ur — dr/dt. 
Equation (13) implies that the specific angular momentum {ru^f, + flr"^) has a constant 
value which, after recalling the initial condition u^po = 0, is readily established to be Qr^. 
This quantity can then be used to evaluate the azimuthal velocity component of the ring, 
yielding the result 

= - nr. (14) 

Similar considerations applied to the r-component of equation (5) lead to the 
equilibrium condition 



GM(r) _ 



y2 J^2 

+ ^ + 2Qu^-- -. ^ = 0, (15) 

r ^ 47rr {pe + p) 
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where we have set — {dur/dt) — 0. Equation (15) expresses the fact that the radial 
equihbrium is characterized by a balance between the outward-directed buoyancy and 
centrifugal forces and the inward-directed Coriolis and magnetic tension forces. Making the 
substitution (r/ro) = 1 + (Ar/ro) and retaining only the lowest order terms in (Ar/ro) 
(<C 1), we obtain 

Ar 1 (GM, 2 1 ( GM, 2peo\ . . 

ro ~8(Qro)H/3oro ""^V 8/3oWHro Peo J ' ^ ^ 

where equation (14) and the approximation pe ~ P [1 + (l//^o)] have been used in simplifying 
equation (15). 

As is evident from equation (16), for a given stellar model, the equilibrium position of 
an equatorial flux ring depends only on the value of the parameter /5o; in the case of the 
present model, (Ar/ro) ~ 6 x 10""^ {(3q/1Q^)~^. This dependence has been confirmed by 
examination of solutions corresponding to a variety of input parameter values, including 
the particular example shown in panels (B), (D), and (F) of Figure 2, and in panels (B) 
and (D) of Figure 3. For this solution, /3o = 10^ as in the case of the ring considered 
above, but oo = 10~^ heo- Because this initial cross-sectional radius is a factor of 10 larger 
than the previous value, the radiative heating rate q and the acceleration produced by the 
aerodynamic drag force are smaller by factors of 100 and 10, respectively (see equations 
[5] and [7]). As can be seen in the relevant portions of Figures 2 and 3, a consequence of 
these reductions is that the ring is subject to vigorous buoyancy oscillations in the radial 
direction. Oscillations of a similar kind in toroidal flux tubes inside non-rotating stars 
have been studied by Spruit & van Ballegooijen (1982), and inside rotating stars by van 
Ballegooijen (1983), Moreno- Insertis, Schiissler, & Ferriz-Mas (1992), and Ferriz-Mas & 
Schiissler (1993). In the present example, the enhanced tendency toward oscillatory motion 
is directly attributable to modiflcations in the buoyant acceleration of the ring; these 
changes arc themselves a result of the altered thermodynamic state of the ring material, 
produced by the diminished rate of radiative heating. The build-up in the magnitude of 
the Coriolis force, together with the action of the drag force, cause the amplitude of these 
oscillations to decrease over time, leaving the ring in the same equilibrium conflguration 
that obtained in the case of the oo = 10~^ /leo ring studied above. 

These results suggest that the equilibrium of a flux ring located in the equatorial plane 
is stable with respect to small amplitude, radial displacements from the position given by 
equation (16). This conjecture has been verified by a series of numerical experiments in 
which an additional, radially directed force of specified amplitude is applied to the ring for 
a brief time interval after it has assumed its equilibrium position. In all cases, the ring is 
initially displaced in radius, but is quickly restored to its original location following cessation 
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of the imposed forcing. Such behavior is not observed in response to the apphcation of 
latitudinal forcing of the same type. In all of these small displacement out of the 

equatorial plane is accompanied by the development of a ^-component of the magnetic 
tension force that accelerates the ring toward the pole. This is similar to the dynamical 
evolution of the so-called polcward-slip instability of equatorial flux rings in the radiative 
layers beneath the solar convection zone (see, e.g., Moreno- Insertis, Schiissler, & Fcrriz-Mas 
1992, and references therein). For the conditions considered in the present paper, the 
perturbed force balance is such that a ring initially contained within the equatorial plane 
moves over time to higher latitudes and larger radii. Once out of the equatorial plane, 
the dynamics of such a ring is found to be identical to that of one with < 7r/2, and is 
discussed in §3.2 below. 

3.2. Flux Rings Outside of the Equatorial Plane 

We now consider the motion of toroidal magnetic flux tubes having initial positions 
To = Rc with < ^0 < 7r/2. We focus on the properties of rings with Oq — n/S (i.e., initial 
latitude Aq = 30°) since the behavior in this case is representative of that displayed by 

all rings having starting locations along the core-envelope interface, out of the equatorial 
plane. In Figure 4, we show the trajectory that a ring with (3q = 10^ and = 10~^ /igo 
follows during the time interval between its release and t = 10^ Ia ~ 10.85 years. As can 
be seen in panel (A), the edge of the ring traces a path that is nearly parallel to the stellar 
rotation axis over this period, extending from tq = 0.232 and Aq = 30° to r = 0.64 R^, 
and A = 72°. In panel (B), we have magnifled the horizontal scale in order to make 
the meridional motion of the ring more apparent. There it can be seen that for a time 
t !v 5 tA after the start of the calculation, the distance r sin 9 of the ring from the stellar 
rotation axis decreases as the ring moves toward higher latitudes along the periphery of the 
convective core. At later times, this moment arm length approaches its initial size, while 
the height r cos 9 of the ring above the equatorial plane steadily increases. 

Additional information pertaining to the evolving dynamical and thermodynamical 
properties of the flux ring under consideration is presented in Figure 5. Inspection of 
panels (A) and (B) reveals that after a brief, initial period of dynamical readjustment, the 
ring moves continuously in the direction of higher latitudes (i.e., the — e^-dircction) and 
larger radii. This is unlike the behavior of rings with Ao = 0° which, in the absence of 
suitable perturbations, do not depart from the equatorial plane as they evolve toward a 
final equilibrium state with Ur = ug = (see §3.1). In the present case, although its upward 
progress slows considerably at later times, the ring never attains an equilibrium in which 
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the forces acting in the e^- and e^-directions balance to produce a state of rest. 

Details concerning the radial and meridional dynamics of the flux ring are provided 
in panels (C) and (D) of Figure 5. When the ring is released at t — 0, it experiences 
a net radial acceleration that is > 0, resulting from the fact that the outward, buoyant 
force is larger than the inward, radial component of the tension force. In the e^i-direction, 

the tension force is initially unopposed by any other force components, leading to a net 
meridional acceleration that is < 0. The ring is therefore pulled toward the pole, and the 
horizontal distance r sin 6 between it and the rotation axis decreases (see Figure 4). This 
change in position has an important consequence for the subsequent motion of the ring. In 
order to conserve angular momentum, the azimuthal velocity increases as the moment 
arm of the ring material becomes smaller (see panel [B]). Associated with this spin- up are 
substantial increases in the magnitudes of the r- and ^-components of the Coriolis force, 
both of which exert considerable influence over the dynamical evolution of the ring. 

Using a procedure similar to that adopted in the case of an equatorial flux ring, the 
azimuthal component of the equation of motion (5) for a ring at higher latitudes can, after 
some manipulation, be written as 



Note that in this case, because the initial radial and meridional motion of the ring 
corresponds to decreasing r sin ^, the sense of the azimuthal motion is such that > 0. 
Hence, the accelerations produced by the Coriolis force in the e^- and e^-directions, 
2 Q Uff, sin 6 and 2 Q cos 9, respectively, are both > 0. As a result, the tension-induced, 
poleward movement of the ring at early times begins to slow when the sum of the 
6'-components of the Coriolis and drag forces becomes large enough to change the sign of the 
net meridional acceleration, making it > 0. Likewise, the acceleration supplied by buoyancy 
and the radial Coriolis force component causes Ur to grow until Ur — —uq cot 9, at which 
point contraction of the ring toward the rotation axis ceases; at later times, Ur > —uq cot 9, 
and the ring expands. The resulting increase in r sin 9 leads to a gradual reduction in u^f,, 
and leaves the ring in a state of near balance between the Coriolis and tension forces in 
both the e^- and e^z-directions. 




(17) 



from which it immediately follows that 




(18) 



The lack of equilibria for flux rings that originate at latitudes Aq > 0° can be 
understood by examining the r- and ^-components of the equation of motion (5) . Assuming 
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that a stationary state with Ur — ug — exists, the r-component can be written in a form 
reminiscent of the equihbrium condition (15) for flux rings in the equatorial plane, 



Pe-P 



GM{r) 2 . 2n 
-O r sm 



\Pe + Pj 

while the ^-component becomes 



+ ^ + 2nu4,sm9- = 0, (19) 

r ^ Anr {pe + p) 



_ f P^\ ^2 ^ 0^!4 ^2 n Ursine + ^ = 0. (20) 

\Pe + pJ r ^ 47rr(pe + p) 

Note that were it not for the appearance of the gravitational acceleration in the radial 
component, equations (19) and (20) would be identical. This correspondence implies that 
the conditions for force balance in the e^- and e^-directions can be simultaneously satisfied 
only if p = Pe- However, as can be seen in panel (E) of Figure 5, p < pe throughout the time 
interval covered by the calculation, indicating that a stationary equilibrium is not possible 
in this particular case. 

In general, all magnetic flux rings with Aq > 0° remain slightly cooler and less dense 
than the surrounding material during most of their ascent, a consequence (in part) of the 
total pressure balance condition given by equation (3). On the basis of equations (19) and 
(20) then, stationary equilibria are ruled out for toroidal flux tubes of this kind. Because 
T < Tg (see panel [E] of Figure 5), a given tube is heated by the radiation that diffuses 
into it from the external medium. In response to this heating and in order to comply with 
the pressure equilibrium condition, the tube expands, becoming less dense as it does so. 
Buoyancy then carries the tube outward to a somewhat larger radius. Hence, apart from 
a brief period at the start of the motion when the buoyant acceleration is established by 
the initial conditions, the rate of rise of the tube is controlled by the rate at which it is 
radiatively heated. This behavior is evident in panel (F) of Figure 5, wherein it can be 
seen that for much of the time interval depicted, the heating and rise times vary in concert. 
Such a situation is in contrast to the dynamical and thermodynamical evolution observed 
previously for flux rings in the equatorial plane. The stationary equilibrium of an equatorial 
flux ring is characterized by T = Tg, a vanishing radiative heating rate, and a balance 
between buoyancy and the sum of the Coriolis and tension forces (see Figures [2] and [3] ) . 

The considerations of the preceding paragraph can be used to derive an approximate 
expression for the radial rise speed of a buoyant flux ring. Interpreting the derivative on the 
left-hand side of equation (7) as (dS/dt) = Ur{dS/dr), the sought-after component of the 
tube velocity is 
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where q and S are defined in equations (7) and (8). To evaluate (d-S/dr), we assume that 
since /3 ^ 1, p ~ Pe and p ~ then, 



5 fa In 



,P0 \PeOj 

where we have made use of the polytrope relation in the form (pe/Peo) = {Pe/PeoY^"- Thus, 

dS /tt — 7\ d In pe fj — a 



dr \ a J dr \ a J he^ 



(23) 



and 



Ur^(-^]qK, (24) 

where he = Pe/ (Ped) is the pressure scale height. In Figure 6, we compare the approximation 
given in equation (24) with the radial velocity component obtained by numerical solution of 
the full equation of motion for j3o = 10^, = 10"'^ heo, and Aq = 30°. At times later than 
t Iai the two results are indistinguishable, indicating that the rate of ascent is indeed 
regulated by the radiative heating of the ring. In fact, using the definitions of the heating 
and rise time scales, equation (24) implies that (x tn, in accordance with the behavior 
seen in panel (F) of Figure 5. 

Figures 7 and 8 contain a summary of results from the simulated rise of a toroidal 
magnetic flux tube released at latitude Aq = 30° with /5o = 10'^ and ao = 10^^ hgQ. As was 
discussed in §3.1, the larger cross-sectional radius of this tube implies that the deceleration 
due to aerodynamic drag (oc a"^) and the radiative heating rate (oc a"^) are both reduced 
from the magnitudes these quantities have when ao = 10~^ /igo- Similar to the equatorial 
flux rings that were the subject of that section, in the present case, these reductions are 
responsible for engendering oscillatory motion of the tube during the early stages of its 
ascent toward the surface. Such behavior can be seen in Figure 7, in which (as in Figure 
4) two views of the path taken by the rising tube are given. In panel (A), it is apparent 
that the larger tube is less buoyant; during the time interval At = 10^ between the 
first and last points on the trajectory, the radial position of the ring only increases from 
To = 0.232 i?* to r = 0.275 i?*. On the expanded scale of panel (B), it is seen that for a time 
fs 0.25 tA following the start of the calculation, the ring moves steadily toward the pole, 
and thereafter executes latitudinal oscillations that decrease in amplitude as the distance of 
the ring from the equatorial plane increases. 

Inspection of the relevant panels of Figure 8 indicates that as in the case of the 
oo = 10~^ heo tube considered earlier in this section, the initial poleward motion is driven 
by the magnetic tension force. However, in the present case, the tube acquires a higher 
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meridional velocity and approaches closer to the rotation axis during this movement, a 
consequence of the fact that the deceleration due to drag is smaller. Also, the change 
in volume arising from the decrease in r sin 9^ together with the diminshed efficiency of 
radiative heating, are sufficient to ultimately make Ap < 0, and to thereby change the sign 
of the buoyant acceleration (see panels [C] and [E]). Note that the the spin- up associated 
with the overall contraction of the ring causes both components of the Coriolis acceleration 
to increase substantiaUy; the oscillatory behavior seen in Figure 7 results from the interplay 
between this acceleration and that produced by the magnetic tension force. 



A primary focus of the present investigation is the buoyant transport of magnetic flux 
from the assumed site of its generation in the convective core to the stellar surface. Of 
particular interest is an estimate of the time required for this transport to take place. Prom 
the definition of Ur given in equation (4), it follows that 



is the time it takes a flux ring to travel outward from tq to any radius r > Tq. 

We have computed tr{r) for toroidal flux tubes having the properties listed in Table 
1. For each case, the radial component of the tube velocity was determined by numerical 
integration of the full equation of motion (5) for a time interval < t < ti, where the 
value of ti for a given solution was chosen to be between 10^ Ia and 3 x 10^ Ia, depending 
upon the magnitude of for that solution. For t > ti, Ur was evaluated using an analytic 
approximation to equation (24), derived in the following way. According to equation 
(24), Ur o: q he at times late enough that any transient behavior stemming from the 
adjustment of the tube to the imposed initial conditions has disappeared. The heating rate 
q depends directly on the temperature difference AT between the tube and its surroundings; 
using the pressure equilibrium condition (3), it can be shown that (AT/Te) S^/(87rpe), 
Pe- Inserting this result in the deflnition of q (see Appendix A), we flnd that 
the dependence of Ur on the physical properties of the tube and the external medium is 
given by (see also Gurm & Wentzel 1967; Parker 1975) 



To further simplify this expression, note that as evidenced by Figures 4 and 7, much of 
the ascent toward the surface takes place with r sin ^ constant (i.e., Ur ~ —ug cot 9). 



3.3. 



Rise Times 




(25) 




(26) 
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Table 1. Properties of Magnetic Flux Ring Solutions 



Number 


/3o 




Cd 


So 


i^(0.96 R^) 










(G) 


(years) 


1 


10^ 


lQ-4.0 


1 


5.77 X 10^ 


1.3010 X 10^ 


2 


10^ 


10-3.0 


1 


5.77 X 10^ 


1.3225 X 10^ 


3 


10^ 


10-4.0 





5.77 X 10^ 


1.3009 X 10^ 


4a 


10^ 


10-4.0 


1 


5.77 X 10^ 


1.3020 X 10^ 


5 


10^ 


10-4.0 


1 


1.82 X lO*' 


1.3011 X 10^ 


6 


10^ 


10-3.5 


1 


1.82 X 10^ 


1.3024 X 10^ 


7 


10^ 


10-3.0 


1 


1.82 X 10^ 


1.3037 X 10^ 


8 


10'' 


10-4.0 


1 


5.77 X 10-'' 


1.3011 X 10^ 


9 


10^ 


10-3.0 





5.77 X 10^ 


1.3013 X 10^ 


10 


10^ 


10-4.0 


1 


1.82 X 10^ 


1.3011 X 10^ 



Note.— All solutions have Aq = 30°. 
'^Virtual inertia is omitted from solution 4. 
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Therefore, from equations (1) and (2),p(f=ipe) oca ^ and B oc a ^ oc Pe, so that Ur varies 
as 

Ur oc (27) 

g 

where we have approximated p ex T ^ Tg. To iUustrate the vahdity of equation (27), 
we have evaluated the constant of proportionahty using the properties of solution 1 (see 
Table 1) at the time t — 50 tA- The rise speed so-derived is shown as a function of time for 
t > 50 tA hy the dotted line in Figure 6. 

The results of our rise time computations are depicted in graphical form in Figure 
9. There we show the time tr{r) required to reach a given radius r within the interval 
0.232 < r < 0.960 i?*, for each of the solutions listed in Table 1. It is apparent that 
rings with smaller values of /3o and oq generally traverse the radiative envelope in less time 
than do rings that are less strongly magnetized and have larger cross-sectional radii (see 
also Moss 1989). Based on the discussion of the preceding subsection, rings having smaller 
values of (3o and Oq are initially more buoyant, and are more likely to remain so because of 
the comparatively shorter time it takes to heat them by inward diffusion of radiation. With 
the exception of the rings corresponding to solutions 7, 9, and 10, the time needed to arrive 
at the stellar surface is less than the estimated main sequence lifetime of the 9 Mq model, 
tus ~ 10^° {M^/Mq)[Lq/ L^) 2.4 X 10^ years. However, note that even the more slowly 
rising tubes of solutions 7, 9, and 10 attain radii in the range (0.92-0.95) within a time 
< tMS- In fact, all of the solutions shown in Figure 9 spend more time crossing the outer 
10% of the stellar radius than they do traveling from r — ro to r — 0.9 i?*. From equation 
(27), it can be seen that since Tg decreases and increases as r approaches i?* in the outer 
envelope, Ur becomes small and the rise time increases accordingly. 

Examination of Table 1 and Figure 9 also reveals that the omission of virtual inertia 
and aerodynamic drag from the equation of motion has little impact on ring dynamics and 
rise times (see solutions 1, 3, and 4). This is because the flux rings considered herein do 
not experience significant, impulsive accelerations, and never achieve large rise speeds for 
prolonged periods of time. Instead, most of the ascent takes place at the slow, quasi-steady 
rate given by equation (24) (or the approximation [25]). This behavior accounts for another 
notable characteristic of the results shown in Figure 9, namely, that the rise time profiles 
of several solutions are nearly identical to one another. From equation (26), it is readily 
seen that the dependence of the rise speed on the intrinsic properties of a given fiux ring is 
Ur OC (a^ P)"^- Consideration of the entries in Table 1 indicates that, as expected, solutions 
having a common profile of tr{r) in Figure 9 are characterized by the same value of the 
parameter combination {al /9o)~^- 

Finally, we note that with the preceding developments, it is possible to derive an 
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approximate expression for the rise time tr{r). To do this, we first use the equation of 
hydrostatic equihbrium in the form dpe/dr = —{g Pe)/c% where Cg is the external sound 
speed, to convert the integration variable in equation (25) from r to Pe] this procedure 
yields 



tr{r) = - j dpe 



c: 



PeO 



,Pe g 




(a-2) 



^) ( — )■ (28) 



PeO 



K'eO 



To obtain the second, approximate equality in equation (26), we have replaced Ur by the 
scaling relation (27) and used the fact that Cg = c^q (pe/Peo)'-" where c^q = (a Peo/peo)- 
The integrand in equation (26) can be simplified by assuming that the opacity follows 
Kramer's law, PeT~^'^. This assumption, together with the polytrope relation 

Te oc p^"^^-* (see Appendix B), allows us to express the integrand as a function of pe only; 
performing the integration we obtain 



tr{r) 



11 a - 13 



-eO 



.go UrO, 




(lla-13) 
2 



- 1 



(29) 



If Uro is determined by evaluating equation (26) at Tq, the rise time estimate (27) becomes 

0.9349 



Ur) = 3.0835 x 10' 



ao 



eO 




- 1 



years. 



(30) 



where the interior model of Appendix B has been used to fix the values of all physical 
quantities at the core-envelope interface. As an example, the dashed line in Figure 9 
represents the rise time profile obtained from equation (30) for [(ao/^eo)^ Po] = 10~^ 



4. Conclusions and Discussion 

The foregoing examination of the physics of magnetic flux rings in the radiative envelope 
of a uniformly rotating hot star has revealed a range of possible dynamical behaviors, 
and has enabled us to estimate the efficiency of buoyancy as a means of transporting 
dynamo-generated fields from the core to the surface. Rings located in the equatorial plane 
can attain a stationary equilibrium in which Ur = 0, 7^ 0, and the combined Coriolis and 
tension forces balance buoyancy. This state is stable against infinitesimal displacements 
of the ring in the radial direction, but unstable when the perturbations are meridionally 
directed, causing the ring position to shift from latitude 0°. Apart from transient episodes 
of oscillatory behavior at the outset of their motion, rings released from latitudes other 
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than 0° move to larger radii along paths that are very nearly parallel to the stellar rotation 
axis. During most of the ascent, the rate of rise of such a ring is controlled by the rate 
at which it is radiatively heated, with a near balance prevailing between the accelerations 
produced by the Coriolis and magnetic tension forces in the r- and 6'-directions. 

Strongly magnetized rings with smaller cross-sectional radii experience larger initial 

buoyant accelerations and have shorter radiative heating time scales. As a result, they 
attain higher rise speeds and require less time to traverse the envelope than do rings with 
weaker fields and/or larger cross-scctional radii (see also Moss 1989). A majority of the ring 
models enumerated in Table 1 reach the stellar surface in a time that is less than the main 
sequence lifetime of the 9 Mq star, the only exceptions being solutions with large values of 
/5o and/or oq. If the field strengths and flux tube sizes considered herein are realistic, this 
result suggests that the field produced by a core dynamo could manifest itself in the surface 
layers of a hot star at a relatively early stage of main sequence evolution. In light of the 
discussion of §3, the strength of the field in a buoyant flux ring that reaches the photosphere 
of the 9 Mq stellar model is estimated to be fa Bq (pe*/Peo) ~ 9.44 (/?o/10^)~^/^ G. 

The preceding conclusions are drawn from results obtained using a simple model for 
a thin, isolated, untwisted flux ring immersed in a non-evolving, rigidly rotating stellar 
interior. Given the exploratory character of this investigation, a more detailed treatment 
of flux tube structure and dynamics is neither warranted nor practical. However, several 
potentially important effects that have been neglected in the present analysis can be 
included within the context of the basic model of §2. Among these are large-scale internal 
circulatory flows, the development (over time) of gradients in composition and angular 
velocity, the influence of radiation pressure on tubes inside more massive stars, mass loss, 
and the consequences of magnetic interference with radiative energy transport (e.g., thermal 
shadows and small-scale flows; Parker 1984). 

Of the effects noted above, rotationally driven, meridional circulation may have the 
most significant impact on the results described in §3. As noted in that section, each of the 
ring models listed in Table 1 spends considerably more time traveling through the thin shell 
0.90 R* <r < 0.96 it!* than in ascending from the starting radius to r = 0.90 (see Figure 
9). In solution 1, for example, the former displacement requires approximately 1.29 x 10^ 
years, while the latter occurs in only 1.38 x 10^ years, a difference in travel time of nearly 
a factor of 100. Wc point out that the inclusion of meridional circulation may decrease 
the time needed to transport magnetic flux through the layers just below the photosphere. 
Estimates of the circulation speed inside upper main sequence stars suggest that while such 
flows proceed quite slowly at great depth, they can become fast in the lower-density region 
near the surface (see, e.g., Kippenhahn & Weigert 1994). Hence, flux tube advection by a 
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circulatory flow may act to supplement buoyant transport in the outermost portion of the 
stellar interior, thereby shortening the time between the production of fields in the core and 
their emergence in the atmosphere. The combined effects of buoyancy and advection will 
be investigated in the next paper of this series. 

We are grateful to Yuhong Fan, M. Maheswaran and Paul Charbonneau for helpful 
conversations, and to Steven Jackson for providing the interior model used in the 
calculations, and to Wendy Mukluk for the rising flux tube figure. We thank the referee for 
a useful question regarding tube heating and for suggesting additional references. JPC is 
also grateful for a visiting scientist position at HAO. 
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A. Derivation of the Ring Energy Equation 

To derive the energy equation for a thin, axisymmetric magnetic flux ring, we begin 
with the flrst law of thermodynamics, written in the form, 

dF+^dlipj" dT- 

where S = (p/p)/(7 — 1) is the specific internal energy of the gas, and (dQ/dt) is the rate 
at which heat is input, per unit mass of ring material. Let Fg be the radiative fiux at the 
surface of the ring; in the case of interest to us here, T < Te, and this flux is directed into 
the ring from the external medium. We assume that the ring is sufficiently opaque that 
all of the flux incident upon it is thermalized within it. For a ring of mass M, the speciflc 
radiative heating rate is then 

dQ_AF^_2F^ 

dt~ M - pa' ^^^^ 
where A — 27ra • 27rr sin 9 is the surface area of the ring. If VTg is the local temperature 
gradient in the surrounding stellar material and n is a unit vector normal to the ring 
surface, the diffusive flux of radiation into the ring is approximately 

16 (7 , _ , 16 (T T3 AT , , , 

Fe^-^ (n-VT,)^- , (A3) 

•J f^e Pe "J Pe ^ 

where AT = (Tg — T), and we have estimated |VTe| ~ {AT /a). 

Exphcit calculation of the derivatives appearing on the left-hand side of equation (Al) 
leads to 

d£ d /1\ p dS 

d^^^dt[-p)=J^pd^' 

where the quantity 

has been referenced to the initial thermodynamic conditions in the ring. Substitution of the 
results given in equations (A2), (A3), and (A4) into equation (Al) yields the ring energy 
equation in the form used herein, 

d^_32 {^-l)aT,^ AT^^ 

dt 3 Ke Pe P 0? Te 

The time scale tn for radiative heating of the ring can be expressed in terms of the rate q 
as tH = [S/{dQ/dt)] = q-\ 
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B. Approximate Stellar Structure 



The distributions of p^, Pe, and Tg used in the calculations described in §3 are 
derived from a model for a spherical, non-rotating star with mass = 9 Mq 
and radius = 3.678 Rq. The star is chemically homogeneous with composition 
X — 0.7131, Z — 0.0169, and contains a central, convective core of mass Mc — 0.308 
and radius Rc — 0.232 R^. The ambient physical conditions at the initial radial position 
ro (= Rc) of the flux ring are Peo = 1.3250 x 10^^ dyne cm~^, peo = 4.7239 g cm~^, and 
Teo = 1.9976 X 10^ K. In order to enhance the precision of computations involving ring 
properties that differ from those of the external medium by very small amounts, it is 
helpful to use an analytic approximation to the background stellar model. We discuss the 
derivation of one such representation in this appendix. 

Unlike the convective envelope of the Sun, the radiative envelope of the adopted 
stellar model contains a significant fraction of the mass of the star, 1 — (Afc/A/*) ~ 0.70. 
Consequently, the gravitational acceleration g{r) in the inner portion of the envelope does 
not vary simply as r~^, but instead depends on the detailed distribution of mass. We 
approximate g{r) within the region between the convective core and the stellar surface by 
using the following piecewise continuous fit to the computed interior model: 

2" 



g{z) = gi exp 



92 exp 




Zo< z < Zi, 

Zl < Z < Z2, 
Z2< Z <l, 



:bii 



where z = r/ R^ and 

9i 



1.80 X 10^ cm s-^ 



1.06 X 10° cm s-^ g2 = 8.90 x 10* cm 8"^ g^ 

Ci = 0.27, Ai = 0.31, 6 = 0.40, A2 = 0.34, 

Zo = 0.232, Zl = 0.422, Z2 = 0.641. 

We assume that Pe and Pe satisfy a polytrope relation of the form p^ oc where the index 
a has a constant value. This assumption, together with the approximation for g given in 
equation (-B1), enables us to integrate the equation of hydrostatic equilibrium in order 
to obtain the structural properties of the stellar radiative interior. The results of this 
procedure are: 



PeO 




6 



Ai 



'zq - 6 



(c-l) 



:b2) 
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in zq < z < zi, where ^{x) 



(2/0F)/die 





is the error function; 



Pel 



1 - 



a 



a 



Pel 92 R* 
. Pel . 



exp 



'6 - zi 



— exp 



(B3) 



in ^1 < ;2 < Z2, where pei = Pe{zi) and Pei = Pe(^i); and, 
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in Z2 < z < 1, where pe2 = Pe{z2) and = Pe{z2)- By virtue of the polytrope relation, the 
pressure and temperature in each of the three intervals are given by 
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where i = 0, 1, 2. 



A polytrope with a = 1.3518 yields a reasonable fit to the properties of the computed 
stellar interior model. Note that a polytropic representation of this kind is not completely 
consistent with the detailed model, a consequence of the fact that the simple assumed 
relation between pe and pe provides only an average description of the entire stellar 
envelope. A more accurate fit can be obtained if the index a is allowed to have a distinct 
value in each of the three intervals delineated above (e.g., ai — 1.400, a2 — 1.330, and 
as = 1.277). However, the single index description is adequate for the purpose of the 
present computations. 
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Fig. 1. — An illustration of a flux tube rising through the radiative envelope of a massive 
star. The tube is assumed to be a geometrically thin circular ring with its symmetry axis 
along the rotational axis of the star. The field in the toroidal tube is assumed to be generated 
by a dynamo process acting at the interface of the convective core and radiative envelope 
zones. The tube is assumed to maintain the enclosed field and the initial mass of the tube. 
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Fig. 2. — The time evolution of selected properties of magnetic flux rings with /3q = 10 
located in the stellar equatorial plane {6 = 7^/2, A = 0°). The panels on the left-hand side of 
the Figure contain results for a ring with initial cross-sectional radius ao = 10^'' heo, while 
those on the right-hand side pertain to a ring with ao = 10~^ /igo- Panels (A) and (B) show 
the normalized radial displacement (Ar/ro) of the ring from its initial position tq. Panels 
(C) and (D) show the velocity components Ur (solid line), uq (dotted line), and it^ (dashed 
hne) in units of the initial Alfven speed uao in the ring. Panels (E) and (F) show the radial 
accelerations (in units of u\q/R^) produced by the buoyant (dotted line), inertial/Coriolis 
(dashed line), drag (dashed-dotted line), and magnetic tension (dashed-double dotted line) 
forces. The net radial acceleration is represented by the solid line. 
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Fig. 3. — The thermodynamic evolution of the magnetic flux rings depicted in Figure 2. 
On the left-hand side of the Figure, f3o = 10^, ao = 10~^ heo, while on the right-hand side, 
= IC^, = 10^^ heo- Panels (A) and (B) show the normalized pressure (solid line), 
density (dotted line), and temperature (dashed line) deficits (Ap, Ap, and AT, respectively) 
of the ring material relative to the external medium. Panels (C) and (D) show the entropy- 
like quantity S (sohd line), the radiative heating rate q (dotted line), the heating time-scale 
tn (dashed line), and the local rise time tji (dashed-dotted fine). 
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Fig. 4. — The trajectory of a buoyant magnetic flux ring in the radiative envelope of a 9 
M© star, as described in §3.2 of the text. Panel (A) shows the path followed by a ring with 

— 10^ and flo = 10~^ heo, starting from a position at latitude Aq = 30° (^o = vr/3) on the 
core-envelope interface. The inner and outer dotted lines denote the core and photospheric 
radii, respectively, and the duration of the time interval between the first and last points on 
the path is 10^ t^. The same trajectory is also shown in panel (B), but with an expanded 
horizontal scale to accentuate changes in the distance of the ring from the rotation axis of 
the star. Beginning in the lower right-hand portion of the panel, the dots along the path 
correspond to the times (t/tA) = 10-^ IQ-^, 10-^ 1, 10, 10^, and 10=^. 
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Fig. 5. — The dynamical and thermodynamical evolution of a magnetic flux ring with 
Po — 10^, Oo = 10~^ heo, and Aq = 30°. All quantities are shown as functions of time, 
measured in units of the Alfven time tA- Panel (A) shows the normahzed radial coordinate 
{r/R^) (solid line) and latitude (A/90°) (dotted line) of the ring. Panel (B) shows the velocity 
components of the toroidal flux tube, using the same format as in Figure 2. Panels (C) and 
(D) show, respectively, the accelerations produced by the r— and ^—components of the 
forces acting on the ring. Individual forces are identified according to the line types used 
in the acceleration panels of Figure 2. Panel (E) shows the contrasts in pressure, density 
and temperature between the tube and its surroundings, while panel (F) shows S and the 
radiative heating rate, together with the local heating and rise times. The formats adopted 
for these panels are identical to those used in Figure 3. 
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Fig. 6. — The time evolution of the radial velocity component for a magnetic flux ring 
with /3o = 10^, oo = 10~^ heo and Aq = 30°. The solid hne corresponds to the radial 
rise speed derived from numerical integration of the equation of motion (5), while the 
dashed line represents the approximate rise speed given by equation (24). The dotted 
line for i > 50 represents the rise speed given by equation (27), with the constant of 
proportionality evaluated using the properties of the numerical solution at f = 50 i^. 



-32- 




Fig. 7. — The trajectory of a magnetic flux ring with /?o = 10^, qq — 10~^ hgo, and Aq = 30°, 
as in Figure 3. The dots on the path in the lower and upper portions of panel (B) indicate 
the positions of the ring at times t = tA and 10 t^, respectively. 
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Fig. 8. — The dynamical and thermodynamical evolution of a magnetic flux ring with 
/?o = 10^, ttQ = 10~^ heo, and Aq = 30°. The format of the figure is the same as that of 
Figure 5, with the exception that panel (E) depicts 10^ times the pressure, density, and 
temperature contrasts, plotted on a linear scale to accommodate the changes in sign that 
Ap undergoes. 
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Fig. 9. — The time required for a magnetic flux ring to travel from tq — 0.232 to a given 
radius r in the interval Tq < r < 0.960 R^. The numbers labelling each curve correspond 
to the models hsted in Table 1. The dashed hne represents the rise time estimate (30), 
evaluated for flux rings having [(ao//ieo)^ /9o] = 10~^. 



